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We study the simplest renormalizable scalar leptoquark models where the standard model is augmented only 
by one additional scalar representation of SU(3) x SU{2) x U(l). The requirement that there be no proton 
decay from renormalizable interactions singles out two such models, one of which exhibits an unusual top 
mass enhancement of the /i — >• e"/ decay rate. We analyze the phenomenology of the model with the unusual 
top mass enhancement of loop level chirality changing charged lepton processes in the light of existing and 
upcoming experiments. Both of the models that do not allow proton decay from renormalizable interactions 
have dimension five operators that, even if suppressed by the Planck scale, can give rise to an unacceptably high 
level of baryon number violation. We discuss symmetries that can forbid these dimension five operators. 
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I. INTRODUCTION 

Currently, the standard model describes most aspects of na- 
ture with remarkable precision. If there is new physics at the 
multi TeV scale (perhaps associated with the hierarchy puz- 
zle), it is reasonable to expect measurable deviations from 
the predictions of the standard model in the flavor sector. 
Amongst the experiments with very high reach in the mass 
scale associated with beyond the standard model physics are 
those that look for flavor violation in the charged lepton sec- 
tor through measurements of the processes, /i — > ej 0] and 
/i — > e conversion [2, 3|, and the search for electric dipole 
moments of the neutron, proton and electron. 

Models with scalar leptoquarks can modify the rates for 
these processes. Simple models of this type have been studied 
previously in the literature, including their classification and 
phenomenology |4-11|. 

Our approach is to first identify the minimal renormaliz- 
able scalar leptoquark models containing one single additional 
representation of SU(3) x 517(2) x U(l) and construct the 
most general renormalizable model without any additional 
constraints on the couplings apart from the usual ones, i.e., 
gauge invariance, Poincare invariance, and locality. Given the 
strong experimental constraints on baryon number violating 
processes like p — > ir°e + , we concentrate only on those scalar 
leptoquark models which don't have baryon number violation 
in perturbation theory. Of course there is baryon number vio- 
lation through non-perturbative quantum effects since it is an 
anomalous symmetry. But this is a very small effect at zero 
temperature. Only two models fulfill this requirement. One 
of those two models gives a top mass enhanced /i — > e*y de- 
cay rate. We perform an analysis of the phenomenology of 
this specific model, including the /i — > e^ decay rate, fj, — > e 
conversion rate, as well as electric dipole moment constraints 
focussing mostly on the regions of parameter space where the 
impact of the top quark mass enhancement is most important. 
For lepton flavor violating processes at higher energies such 
as t — > /j,-f, deep inelastic scattering e + p — >• fi(r) + X, etc., 
the impact on the phenomenology of the top quark mass en- 
hancement of charged lepton chirality flip is less dramatic and 
that is why we focus in this paper on low energy processes 
involving the lightest charged leptons. 



We also consider the effects of dimension five operators that 
can cause baryon number violation. We find that the two mod- 
els without renormalizable baryon number violation can have 
such operators and, even if the operators are suppressed by 
the Planck scale, they may (depending on the values of cou- 
pling constants and masses) give rise to an unacceptable level 
of baryon number violation. We discuss a way to forbid these 
dimension five operators. 



II. MODELS 

A general classification of renormalizable leptoquark mod- 
els can be found in [4] [5). However, in the spirit of our ap- 
proach, in which we are interested in models with no proton 
decay, a more useful list of possible interaction terms between 
the scalar leptoquarks and fermion bilinears is presented in 
fl2l . where those models that have tree-level proton decay are 
highlighted. The relevant models are listed in Table I below. 



leptoquark 
couplings 


diquark 
couplings 


517(3) x 51/(2) x (7(1) 
representation of X 


XQe, XLu 


- 


(3,2,7/6) 


XLd 


- 


(3,2,1/6) 


XQL, Xue 


XQQ, Xud 


(3,1,-1/3) PD 


XQL 


XQQ 


(3,3,-l/3) PD 


Xde 


Xuu 


(3,l,-4/3) PD 



TABLE I: Possible interaction terms between the scalar leptoquarks and 
fermion bilinears along with the corresponding quantum numbers. Represen- 
tations labeled with the subscript "PD" allow for proton decay via tree-level 
scalar exchange. 

The only two models fulfilling our requirement are 

X = (3, 2, 7/6) and X = (3, 2, 1/6) . 

Model I: X = (3,2,7/6). 

The Lagrangian for the scalar leptoquark couplings to the 
fermion bilinears in this model is, 



£ 



\^ R X T eV L - A«gUftQ£ + h.c. , (1) 



where, 



X 





L, 




(2) 



After expanding the SU(2) indices it takes the form, 
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Y a v{) 



h.c. 



(3) 



Note that in this model the left-handed charged lepton fields 
couple to right-handed top quarks, and the right-handed 
charged lepton fields couple to left-handed top quarks. So a 
charged lepton chirality flip can be caused by the top mass at 
one loop. 

Model II: X = (3,2,1/6). 



The corresponding Lagrangian is, 
C = -X^^X T eL{ 



(4) 



where we have used the same notation as in the previous case. 
Expanding the SU(2) indices yields, 



C 



-\ l i^ H {V a e 3 L 



A d u aR\ 



Y a <) 



h.c. 



(5) 



In model II the leptoquark cannot couple to the top quark, so 
there is no m* enhancement in the /j, — > cy decay rate. There 
is also no m& enhancement, and the one-loop effective Hamil- 
tonian for ^ — } cy (after integrating out the massive scalars 
and the heavy quarks) is proportional to the muon mass. For 
this reason, in the remainder of the paper we will focus en- 
tirely on model I. 



III. PHENOMENOLOGY 

In this section we analyze some of the phenomenology of 
model I, i.e., X = (3, 2, 7/6). We concentrate only on those 
constraints which are most restrictive for the model and po- 
tentially most sensitive to the unusual top mass enhancement 
of the charged lepton chirality change, i.e., the ones coming 
from the following processes - muon decay to an electron and 
a photon, muon to electron conversion, and electric dipole mo- 
ment of the electron. 



A. /j, — > ej decay 

The relevant Feynman diagrams for this process are pre- 
sented in Fig. 1 . The uniqueness of model I is that, apart from 
the fact there is no tree-level proton decay, the /i — > ej rate is 
enhanced by the top quark mass. To our knowledge, such an 
enhancement of /j, — > cy was observed previously only in [8 1 
in the context of an SU(2) singlet scalar leptoquark model. 
However, that model suffers from perturbative proton decay 
and the impact of the m t enhancement was not focussed on. 





FIG. 1: Feynman diagrams contributing to the process n — > ey. 



Keeping only the piece enhanced by m t , the sum of ampli- 
tudes corresponding to the diagrams in Fig. 1 (neglecting the 
terms proportional to m e ) is given by, 



iM = 



em t 



16ir 2 my 



f{mymy)k v e^{k) 



+ (\l 1 yC\?ye L (p-k)a^Mp) 



(6) 



where k is the photon four-momentum and e is the photon 
polarization. The function f{m 2 /m 2 / ) is given by, 
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(7) 



and the tilde over the couplings denotes that they are related by 
transformations that take the quarks and leptons to their mass 
eigenstate basis through the following 3x3 matrix transfor- 
mations, 

\ u = U{u,R)^\ u U{e,L) , A e = U(e, JE)+ X e U(u, L) , (8) 

where the right-handed up quarks in the Lagrangian are re- 
lated to the right-handed mass eigenstate up-type quarks by 
the matrix U(u, R), the left-handed up quarks in the La- 
grangian are related to the left-handed mass eigenstate up-type 
quarks by the matrix U(u, L), etc.. 
The fi—} e-f decay rate is, 



T(fM -> ej) 



e 2 X 2 m 2 m^ 
2048 7T 5 mt 



f(mi/m 2 v y , 



where, 



A = 



1 i~ ~ i2 1 i~ ~ ,2 

_ A 13 A 32 + - A 31 A 23 



(9) 



(10) 



Fig. 2 shows the relation between A and the scalar leptoquark 
mass. This dependence was plotted for the /i — > ej branching 
ratio equal to the current upper limit of Br(/i — > ey) ~ 2.4 x 
10~ 12 reported by the MEG experiment, and the prospective 
MEG sensitivity of Br(/i -> ey) ~ 5.0 x 10~ 13 . It shows that 
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FIG. 2: The combination of couplings A from Eq. \\()\ as a function 
of the scalar leptoquark mass for two values of the /i — > ej 
branching ratio relevant for the MEG experiment. 



the experiment will be sensitive to scalar leptoquark masses at 
the thousand TeV scale for small values of the couplings. 

For very small x, f(x) — > f(x) = | log a;. This is a rea- 
sonable approximation in the range of x we are interested in. 
For example, /(lCT 8 )//(lCr 8 ) ~ 1.1. 



B. fi — > e conversion 

The effective Hamiltonian for the \i — > e conversion arises 
from two sources, 



rtoff — n c g + n cS 



(11) 



The first is the dipole transition operator that comes from the 
loop diagrams which are responsible for the /j, — » ej decay, 
given by, 
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Using the following Fierz identities (for spinors), 

{u1LU2r){u 3 H.U4l) = ^(uiL7 M W4L)(w3i?.7M u 2fl) , 
{u1LU2r){u 3 lU4:R) = ■7 ) {uiLUiR){u 3L U 2 R) (13) 

+ ■^{u\L'7 tlU U A R){u 3 R(J l _ l yU2L) , 

we arrive, after integrating out the heavy scalar leptoquarks 
(at tree level), at the second part of the effective Hamiltonian, 
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iff = Y^~> ^ui^uYizLl^VL^UaRl^UaR) + A^A* 2 C S (^)(eRfi L )(u a RU aL ) + j Ct (/ti) (eflO"'"' 'l± L ){u a RO nvU aL ) 



+ X e (K )*{.e-R.l tl ^B){ u aLl^U aL ) + {\ e )* (A u )* C S (lJ,)(e L fl R )(u a LU a R) + -7C T (n){e L (T ,i '' ' (iR)(u>aL<rnvUaR) 



+ ^-^-(XeVcKMY 1 {(XeVcKM) l ) {zrY ^R){d a Ll^d aL ) + 



(14) 



The CKM matrix arises whenever a coupling to the left- 



handed down- type quark appears. In Eq. ( 14 1 the contribu 



tion of the heavy quarks, as well as the contribution of the 
strange quark, are in the ellipses. Since the operators qq and 
qa^q do require renormalization, their matrix elements de- 
velop subtraction point dependence that is cancelled in the 
leading logarithmic approximation by that of the coefficients 
Cs.t- Including strong interaction leading logarithms we get, 



and 
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a s (n) 
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-12/(33-2iV,) 



-|4/(33-2JV„) 



(15) 



(16) 



where N q = 6 is the number of quarks with mass below my- 
In order to match the effective Hamiltonian ( fl4| i to the Hamil- 



tonian at the nucleon level and use this to compute the conver- 
sion rate, we follow the steps outlined in |[T3l[T4ll . 

Our results, taking into account only the contribution from 
T-lgg, are shown in Fig. 3. The current experimental limit is 
Br(/i — > e conversion in Au) < 7.0 x 10~ 13 fl5l . However, 
here we focus on the prospective Mu2e experiment [2 1, which 
has a sensitivity goal of 5 x 10~ 17 . The COMET experiment 
1 3 1 aims for comparable sensitivity in later stages. We use the 
total capture rate for f|Al of w capt uro = 0.7054 x 10 6 s _1 
ifTBI to switch from the fi — » e conversion rate to a branching 
ratio. 

Apart from coupling constant factors, the contribution to 

(a) 

the /i —$■ e conversion amplitude from rL e g is enhanced 

over the contribution to the amplitude from 7i^ s roughly by 
(TO t /TO M )(e 2 /327r 2 )log(m 2 / /m 2 ) ~ 7, for m v in the thou- 
sand TeV range. 

Our results show that in some regions of parameter space 
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FIG. 4: The combination of couplings A from Eq. \\i)\ as a func- 
tion of the scalar leptoquark mass for a branching ratio Br(/i — > 



e conversion in Al) = 10 and four different positive values of 



FIG. 3: The combination of couplings A from Eq. dTOb as a func- 
tion of the scalar leptoquark mass for two values of the Br(/x — > 
e conversion in Al) relevant for the Mu2e experiment. The dashed 

line, corresponding to Br(/i -> ej) = 1CT 14 , is included for refer- the ratio of the couplings r from Eq. (TTJ 
ence. 



the Mu2e experiment will be able to constrain leptoquark cou- 
plings with similar precision to what can be done with an ex- 
periment which is sensitive to a branching ratio for /i — > e^ 
of around 10 -14 . In other regions the Mu2e experiment is 
likely to give a more powerful constraint for such a /j, — » e^ 
branching ratio, for example, when the Yukawa couplings 
are strongly hierarchical and the top quark loop is very sup- 
pressed. 

To show graphically the contributions to the branching ra- 
tio originating from terms in the effective Hamiltonian with 
different structures, we set all the couplings to zero apart from 
Ag 3 , A 23 , A 31 , A 32 , A* 1 , A^ 2 for simplicity, i.e., we leave only 
the couplings relevant for the [i — » ej decay and one of the 
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vector contributions to H 

Note that the heavy quark contributions are suppressed 
by Aqcd/wq, low energy phenomenology suggests that the 
strange quark contribution is small, and furthermore the ten- 
sor contributions are not enhanced by the atomic number of 
the target. 

In addition, we consider only real couplings and define k = 
A* 1 A* 2 . We also assume A;! 3 A 32 = A 31 A 23 = A, so that we 
can plot A as a function of the scalar leptoquark mass my for 
a given value of the ratio, 



11\12 



\}\ 



i(A e 3 A 32 ) 2 + i(A 3 iA 23 )s 



(17) 



Figs. 4-7 show our results for a few values of r = ±1, 
±7, ±25, ±50 and two values of the branching ratio Br(/i — ¥ 
e conversion in Al) = 10~ 16 , 10 -17 . 

fa) 

For r < 1 the branching ratio is dominated by the W cf { 
contribution and in this parameter region all curves look like 
the ones in Fig. 3. For larger values of r, depending on the 
relative sign between the contributions from Hj^ and n cS , 
there are two possibilities. If the interference is constructive, 
the curve moves down with increasing r since a smaller value 
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FIG. 5: Same as Fig. 4, but for negative values of r. 



of the coupling A is required to achieve a given branching ra- 
tio (Figs. 5, 7). In the case of a destructive interference, the 
curves move up until a value of r is reached for which the two 
contributions are the same (Figs. 4, 6). As estimated before, 
this occurs for r sw 7. Increasing r further brings the curves 
back down, since the 'H e J contribution becomes dominant. 

Large values of r are expected if the Yukawa couplings of 
X exhibit a hierarchical pattern like what is observed in the 
quark sector; k changes generations by one unit while the 
product of couplings in A involves changing generations by 
three units. Finally, we note that for all the curves in the plots 
above the Yukawa couplings are well within the perturbative 
regime. 



C. Electron EDM 

Another flavor constraint on the couplings of model I comes 
from the electric dipole moment (EDM) of the electron. As 
mentioned earlier, the fact that X couples directly to both left- 
and right-handed quarks means that at one loop the top quark 



0.12 


Br(//->e conv.in 


Al)= 


=10" 


-17 




/ '- 


0.10 


- 








r= 


=7 / / 


0.08 


- 










rr=\S 


0.06 


- 










S 


0.04 












r=25^^: 


0.02 


^ 










r^o,.: 


00 












; 



2000 



4000 6000 

m v [TeV] 



8000 




A" 




Q,u 



e,L 



FIG. 8: Feynman diagram representing proton decay in model I. 



IV. BARYON NUMBER VIOLATION AND DIMENSION 
FIVE OPERATORS 



FIG. 6: Same as Fig. 4, but for a branching ratio Br(/i 
e conversion in Al) = 10 -17 . 
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FIG. 7: Same as Fig. 5, but for a branching ratio Br(/^ 
e conversion in Al) = 1CP 17 . 



mass can induce the chirality flip necessary to give an electron 
EDM. We find that, 




emt 



13\31i 



\6-K 2 rriy 



Hmi/mWMK'K 



(18) 



The present electron EDM experimental limit [ 17 1 is, 



|4| < 10.5 x l(T 28 ecm. (19) 



We can write the dipole moment in terms of the branching 
ratio, Br(/i — > erf), giving the constraint 



A 



VBr(/x->-e7) < 2.0 x 1(T 7 



(20) 



For example, if model I gave a branching ratio equal to the 
current experimental bound of Br (/i — » ej) < 2.4 x 10~ 12 , 
this would correspond to the constraint on the couplings of 

|lm[Ai 3 A3 1 ]|/A<0.13. 



Tree-level renormalizable interactions are not the only pos- 
sible source of baryon number violation. It might also occur 
through higher-dimensional nonrenormalizable operators. In 
the standard model proton decay is restricted to operators of 
mass dimension six or higher. However, the scalar leptoquark 
models we consider exhibit proton decay through dimension 
five operators. 

Let's first consider model I, in which X = (3, 2, 7/6). Al- 
though it doesn't give proton decay at tree level, one can con- 
struct the following dimension five operator, 



Or 



9 ab d a Ra d b Rp (H^X y 



a/37 



(21) 



The coupling constant matrix g is antisymmetric in flavor 
space. Because of the tree-level leptoquark couplings (see, 
Table I), baryon number violating decay occurs here through 
the process shown in Fig. 8, resulting in n — >• e~K + and 
p — > K + v. Setting the coupling constants to unity, we esti- 
mate the baryon number violating nucleon decay rate caused 
by this operator to be, 
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(22) 



Since the current experimental limit is Tp Xp < 2.7 x 
10~ 66 GeV [18 1, even if the scale of new physics A is equal 
to the Planck mass Mpl when the coupling constants are 
unity, this operator causes too large a proton decay rate for 
m v < 10 000 TeV. 

In the case of model II, where X = (3, 2, 1/6), there are 
two dimension five baryon number violating operators, 



O 



(i) 
ii 



g ab u Ra d b R ^X 1 



a/3-r 



Of} = j9 ab u Ra e R (X^X^)e a ^ 



,(i) 



(23) 



The operator O}/ permits a nucleon decay pattern similar to 



the previous case, e.g., 



and p — > 



Proton 



(2) 

decay through the operator 0\J is much more suppressed. 

In order to prevent proton decay through dimension five op- 
erators, one could introduce a discrete gauge symmetry that 
forbids the baryon number violating nonrenormalizable cou- 
plings. Since B — L is the only anomaly free global sym- 
metry in the standard model, we chose to impose a discrete 



subgroup of B — L. In models I and II the leptoquark has 
B — L = 4/3. The usual Z2, where the nontrivial transfor- 
mation is (— 1) B ~ L , doesn't work, as the operators Oi, O i1 , 

and Ojj are invariant under this transformation. However, we 
find that imposing a Zj, discrete symmetry, with elements that 
are powers of exp[2iri(B — L)/3], forbids these dimension 
five operators and, thus, prevents the proton from decaying 
in this class of models. Note that gauging B — L and spon- 
taneously breaking the symmetry with a charge three scalar 
(at some high scale) leaves this unbroken discrete Z3 gauge 
symmetry. It is not possible to use any discrete subgroup of 
B — L to forbid proton decays in the models from Table I 
which exhibit proton decay at tree level since all the interac- 
tions conserve B — L. 

Finally, we would like to comment on the relation between 
this work and that of [12|, where renormalizable models that 
have additional scalars and have baryon number violation at 
tree level but not proton decay were enumerated and dis- 
cussed. In these models none of the scalars were leptoquarks 
(they could rather be called diquarks or dileptons). However, 
if we permit higher dimension operators, then models 4 and 9 
containing the scalar X = (3, 1, 2/3) (which has renormaliz- 
able diquark couplings), have dimension five leptoquark-type 
couplings, 



O 



hi 



g a \QTH)e b R X 



(24) 



This operator, combined with the renormalizable couplings of 
X to two quarks, gives proton decay with the rate estimated 
in Eq. (22 1. This observation restricts the parameter space of 



models 4 and 9 presented in 1121 to the one in which either the 
color triplet scalar X is very heavy or its Yukawa couplings 
are small. 



V. CONCLUSIONS 

We have investigated the minimal set of renormalizable 
models in which a single scalar leptoquark is added to the 



standard model with the requirement that proton decay not be 
induced in perturbation theory. We have looked in detail at 
one particular model which gives an unusual top quark mass 
enhancement of the branching ratio of p, — > cy. 

For this model, we have compared the /i — > ej branching 
ratio to the fj, —} e conversion rate in light of current con- 
straints and future experiments. We find that the most strin- 
gent constraints on this model could come from the Mu2e 
experiment, and we have shown the potential limits both the 
MEG and Mu2e experiments could place on some of the cou- 
plings of the scalar leptoquark to the Qe and Lu bilinears. 

We have also shown the constraints on this model coming 
from the most recent limits on the electron EDM. Although 
the electron EDM also has the unusual enhancement from the 
top quark mass, the constraints are not so strong when com- 
pared with lepton flavor violating effects. 

We have commented on the existence of nonrenormalizable 
operators in these minimal models which can give an unac- 
ceptably large proton decay rate for ray < 10 000 TeV, as 
well as provided a simple mechanism for avoiding them. 

Since there are only two scalar leptoquark models where at 
the renormalizable level baryon number is automatically con- 
served, it would be interesting to examine a more extensive 
range of phenomena and address, over a wide range of param- 
eter space, how to distinguish experimentally between these 
two models. 
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